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The charge and spin patterns of a quantum dot embedded into a spin-orbit coupled quantum 
wire subject to a magnetic field are investigated. A Luttinger liquid theory is developed, taking 
into account open boundaries and finite magnetic field. In the quasi-helical regime, when spin-orbit 
effects dominate over the Zeeman interaction, peculiar states develop at the Fermi surface of the dot. 
Anomalous Friedel oscillations with twice the expected wavelength develop in the wavefunction of 
collective excitations of such states, accompanied by peculiar spin patterns in their magnetization. 
Both effects are analyzed in detail and shown possible to be probed in transport experiments. 
The stability against electron interactions and magnetic field is investigated. We also discuss how 
signatures of such states survive in the total charge and spin densities. 
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I. INTRODUCTION 

When electrons are confined in a tight one-dimensional 
portion of space, namely in a one-dimensional qua ntu m 
dot, marked oscillations occur in the charge densitj^^^. 
Friedel oseillation^ develop due to the presence of con¬ 
fining barriers which break the translational invariance 
and induce backscattering at the edges. Such oscillations 
exist regardless of the interactions between electrons, and 
give rise to a number of peaks in the charge density pro¬ 
portional teP Nq/D, where Ng and D are the total num¬ 
ber of electrons and the degeneracy of single particle lev¬ 
els respectiveljP^. As an example, for a one-dimensional 
quantum wire of spinful electrons one would expect No/2 
peaks for even Nq. 

On the other hand strong interactions among the parti¬ 
cles lead to the formation of peculiar correlated states, 
dubbed Wigner moleculeJ^^®^, the finite-size counter¬ 
part of Wigner crystalP. Such a molec ule is character¬ 
ized by Nq peaks in the electron densitjPEEl, regardless 
of £>. 

Typically, in all systems with a degeneracy D > 1, Friedel 
and Wigner oscillations have different wavelength and 
can be thus discriminated by looking at the charge den¬ 
sity. On the other hand, in systems with D = 1 - such 
as for instance a spin-polarized one-dimensional electron 
liquid - Friedel and Wigner oscillations have the same 
wavelength and cannot be distinguished at the level of 
single-particle densitj^i^. 

Many of these effects have been investigated in one¬ 
dimensional semiconducting quantum wireJ^^ and car¬ 
bon nanotubeJi^, in which a quantum dot can be created 
by defeetJi^, suitably crafted tunneling barriers or even 
buckling a carbon nanotub^^. 

Recently, novel peculiar one-dimensional systems have 


been created, the helical liquids, occurring for example 
at the edges of topological insulators^^ltlSl or in carbon 
nanotubes subject to an electric fielcP^. In a helical 
liquid, electrons with opposite spin counter-propagate 
due to spin-momentum locking. In the presence of 
time-reversal (TR) symmetry, spin-momentum locking 
protects the chiral propagation of electrons and prevents 
elastic backscattering: indeed, non-magnetic barriers are 
not effective in confining the system. On the other hand, 
magnetic barriers can induce backscattering, leading to 
the formation of spin density waves^^ Two such barriers 
can create a quantum dot into the helical system, in 
which peculiar spin textures and spin ordering occurs 
and can be controlled by means of static or AC magnetic 
perturbations However, such magnetic barriers do 

not give rise to charge oscillations in stark contrast with 
the previous case. 

Spin-orbit coupled quantum wires subject to a magnetic 
fieicp^OH gjgo systems which have been lately in¬ 
vestigated in depth. They exhibit a quasi-helical (Q-H) 
behavior different from the helical liquid discussed above. 
Indeed, the magnetic field breaks the TR symmetry 
mixing left- and right-movers and opens a Zeeman gap 
at zero momentum (see Fig. [^. States in this gap 
display a peculiar character, due to the non-perfect 
spin-momentum locking. Since in these wires TR is 
broken by the magnetic field itself, a quantum dot can 
be created by means of usual non-magnetic barrier^^HHlI] 
In such a dot, charge oscillations will occur and due 
to the existing correlation between chirality and spin, 
peculiar spin textures are expected. 

Such nanowires have been widely investigated in the 
last few years. When proximized with an s-wave su- 
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perconductor, Majorana states occur at the boundaries 
of the wir^^, which can be in principle investigated via 
STM transport!^. Also, an enhancement of the gap in¬ 
duced by electron interactions has been reportecP^, with 
anisotropic spin propertie^^ and spin textures in the 
presence of magnetic impuritie^^. Quasi-helical states 
also occur in wires with hyperfine coupling to the nu¬ 
clear spinJ^. 

Such theoretical studies are based on a Luttinger 
modeP^f^ developed in the limit of vanishing applied 
magnetic held, typically employing periodic boundary 
conditions. 

The task of this paper is to directly investigate the Q-H 
states which develop at hnite magnetic held within the 
gap of a quantum dot and to assess how they affect the 
charge and spin densities. To do so, we consider spin- 
orbit and a non-vanishing magnetic held on equal foot¬ 
ing. Employing open-boundary conditions, states within 
the band gap are considered, for which a linearized spec¬ 
trum and the corresponding wavefunctions for the single¬ 
particle problem are obtained. A Luttinger model with 
open boundaries is then developed, valid when the Fermi 
energy lies within the band gap, also in the presence of 
interactions among the electrons. 

We study in details the charge distribution and the mag¬ 
netization of states near the Fermi surface, introduc¬ 
ing the concept of collective excitation wavefunction and 
collective excitation magnetization. They are the one¬ 
dimensional analogue of the quasiparticle wavefunction 
already introduced in literature and can be probed by 
means of STM transport experiment d^° * ^^ l The stability 
of Q-H states against the intensity of the applied mag¬ 
netic held and electron interactions is assessed. 

We also investigate how the properties of such Q-H states 
rehect on those of the total charge and spin densities. 
Since these quantities also involve states below the Fermi 
surface one needs to go beyond the Luttinger theory. To 
do so, we have employed a Hubbard model solved by 
means of an exact diagonalization procedure in the ab¬ 
sence of interactions and by a variational MPS algorithm 
in the interacting regimd^2I43i 

Our main results are the following. When the spin-orbit 
effects dominate over the magnetic held, the Q-H states 
exhibit peculiar charge oscillations. In stark contrast to 
what one would naively expect for states with D = 1, 
they do not exhibit Nq peaks as for conventional Friedel 
oscillations, but rather Nq/2. We dub this an anomalous 
Friedel oscillation and show that it is intimately con¬ 
nected to the presence of evanescent states which form 
at the edges of the quantum dot within the band gap. In 
the Q-H regime, the length scale of these states becomes 
comparable to the dot size, resulting in the formation of 
anomalous Friedel oscillations. The latter are accompa¬ 
nied by peculiar spin textures: although the magnetiza¬ 
tion of the Q-H states precesses with Nq peaks, strong 
modulations of the magnetization modulus occur near 
the dot edges, resulting in an effective doubling of the 
wavelength. Increasing the magnetic held or the inter¬ 


action strength results into a progressive disappearance 
of the anomalous Friedel oscillations and ensuing spin 
textures with the dot edges being more stable with re¬ 
spect to its center. Signatures of the anomalous Friedel 
oscillations can be detected also in the total charge den¬ 
sity, although the effects in this quantity are much less 
striking. 

The paper is organized as follows. In Section [IT] we intro¬ 
duce the model and construct the Luttinger liquid the¬ 
ory in the presence of open boundary conditions. In Sec- 
tion lHIl we evaluate the collective excitation wavefunction 
and magnetization, describing the anomalous Friedel os¬ 
cillations and the peculiar spin textures. We also discuss 
their stability against the magnetic held and electron in¬ 
teractions. Finally, we study the charge and spin densi¬ 
ties. Section HVl contains the conclusions. 


II. THE MODEL 

A. Single-particle problem 

We consider a quantum dot of length L, with a Rashba 
type spin-orbit interaction r] = —rjUz (rj > 0) and subject 
to an external magnetic held B — Bu^ {B > 0), with 
and Uz axis unit vectors (Fig. [^. Its single-electron 



FIG. 1. Representation of the quantum dot of length L with 
intrinsic spin-orbit interaction, characterized by the vector r], 
subject to an external magnetic held B, parallel to the wire 
axis and orthogonal to the spin-orbit interaction. 


Schrodinger equation (fi, = 1) is 

(~2^ + iWzdx + ^g*p.BBa^ = E^{x). (1) 


Here p-B the Bohr magneton, g* and m* are the ef¬ 
fective gyromagnetic factor and the band mass of the 
electrons respectively, ai are the Pauli matrices and 
'I'(a:) = is the spinor eigenfunction satis¬ 

fying open boundary conditions (OBC): 'I'(O) = '^{L) = 
0. To describe the system the following dimensionless 
parameters are useful 


a = 


FnL’ 


P = 


9*1^bB 
2En ’ 


( 2 ) 
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FIG. 2. Band structure of the quantum dot in the Q-H regime 
with k in units L~^. Parameters: a — 10^, /? = 10^ {5 = 0.1). 


where Eq = {2m*LF')~^. In the rest of the paper all the 
energies will be written in units of Eq. The spectrum 
obtained from Eq. Q consists of the two bands 

e±{k) = (3) 


where k G with a set of discrete wavevec- 

tors determined by the OBC, to be specified later. The 
magnetic field opens a gap A = 2/? at /c = 0. The pa¬ 
rameter 


S = 




( 4 ) 


ide ntifies tw o opposite regimes: the Q-H one for S < 
l/2j ^6 | 27 | ^ dominated by spin-orbit, and the one for 
6 > 1 / 2 , where conversely the external magnetic field 
is prevalent and the system begins to polarize. When 
5 <1/2 the band S-{k) has a maximum at fc = 0. The 
wavevectors inside the gap rarige from k- to fc_|_, with 
k± = L~^ ±2/3 (see Fig. [^. On the other hand, 

when 5 > 1/2 the two minima vanish and the band £_ (k) 
exhibits a single minimum in fc = 0 (not shown); in this 
case, fc_ = 0. Equation ([^ can be formally solved for k 
as a function of e with the result 


fci.3(e) 


1 We E cP- E ^-|- ^ 

ZV 2 


( 5 ) 


where the -I- and — signs refer to ki and fca respectively. 

In this paper we are interested in studying energies 
inside the gap, i.e. |e| < /3, where kiL G K+ and k^L is 
a purely imaginary number. It is then useful to rewrite 
it as fca = ZKs, with kqL G K+. Thus, the most general 
single-particle spinor wavefunction satisfying Eq. Q for 
|£| < /3 has components 


V’T(a^) 




y2L 

( 6 a) 

y/2L 

( 6 b) 


with 


, £ — -|- aLki 

di= ^ Cl, 

(7a) 

£ — — Laki 

d2= ^ 

(7b) 

, £ + -|- iaLn3 

d3= ^ 

(7c) 

£ + — iaLK3 

«4 = ^ C4 , 

(7d) 

obtained from Eq. Q. 

Imposing OBC, from Eqns. Q and Q 

one gets a sys- 


tern of four linear equations that can be written as 

M 0=0, ( 8 ) 

with M a 4 X 4 matrix of coefficients (not reported) and 
c = (ci, C 2 , C 3 , 04 )^. Imposing the condition Det(M’) = 0 
we obtain the secular equation 


a^\JPP — £^ [cos(A:iL) cosh(K 3 L) — 1] = 

(a^£ -I- 2/?^) sin(fciL) sinh(K 3 L), (9) 

which corresponds to an implicit equation for e. 

In the following we will consider the solution of this equa¬ 
tion in the limit K 3 L ^ 1. Indeed, as shown in Ap¬ 
pendix 1 ^ this corresponds to have a large number of 
states in the gap, i.e. — fc_ 3> tt/L, a necessary re¬ 
quirement in order to construct a Luttinger theory for 
states within the gap (see Sec. IIB). Employing this con¬ 


dition and concentrating on the states near the center 
of the gap (|£| <C /3) one obtains analytical expressions 
for the linearized spectrum ^ and for the discretized 
wavevectors and K 3 ^„ (with n G Z), see Appendix [A| 
The main results are: 




■^0 



( 10 ) 



TT 



( 11 ) 


where k^'^ = fci(£ = 0) = aL ^[(\/l -|- 46^ + l)/2]^/^ is 
the reference wavevector. 


oP\/l E 4^2 


( 12 a) 


the corresponding velocity and 


2 

7=1 -arctan(2i5). (12b) 

TT 

We choose as a reference point n « no, with £„j, « 0 the 
closest level to £ = 0 . As discussed in Appendix [A] one 
can always approximate K 3 _„ as 


^3,n 


.(0) 


= K3(£ = 0) = -j 


OL ( — 1 


1/2 


L \ 


(13) 
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The coefficients ci, ...,C4 can be written as a function 
of Cl only via the system in Eq. (|^. For |e| <C /3 and 
» 1 they are 


C2= • 


X L 


\/2(x+l)-l 


^ i Vx+l-v^ 

X Vx-1 


1 

C3 = - 
X 


C4 = (-1) 


v^X - Vx + 1 , 


Cl, 


n+1 


g— 


Cl 


Cl, (14a) 
(14b) 
(14c) 


where y = Vl + 4(5^. Coefficients are then ob¬ 

tained from Eq. 0 and Cl € K is numerically determined 
from the normalization of the spinor. 

The above relations are valid in the whole range of 
5. However, the polarized regime (<5 > 1/2), which can 
also be discussed in the framework of a spinless Luttinger 
liquicP^, is well known and not particularly novel. Thus, 
in this paper we will focus mostly on the Q-H regime, 
5 < 1/2. We emphasize here that even in the regime 
5 ^ 1 / 2 , which is the most interesting one, we continue 
to assume the presence of a sufficiently large number of 
states in the gap with a finite magnetic field. This is 
possible since a and /3 are independent parameters so 
that one can always tune 6 = fi/a^ ^ 1/2 still satisfying 
the condition Kg^^i = ft/a » 1, see Appendix |a[ In 
addition, in the Q-H regime the ratio between the number 
of states in the gap. A/,, compared to the total number 
of states in the Fermi sea. A/), is Ng/Nt oc 5. Thus in the 
Q-H regime the Fermi surface properties of the system 
are well described by our theory, conversely the physics 
of the whole Fermi sea is not necessarily captured. 

For future convenience analytical expression of the co¬ 
efficients up to first order in 5 are quoted here: 


Cl = 1, (15a) 

C 2 = (15b) 

5 , , 

C3 = -1 + * 2 > 

C4 = (-l)"+ie-'=3'”^ + ^6^ , (15d) 



d2 = i + 2S, 


= -(. + ?«), 

f-i 



(16a) 

(16b) 

(16c) 

(16d) 


if^{x) the leading term is the right-moving one, oc 
while in tfiix) is the left-moving one, oc In par¬ 

ticular, at zero-order in <5 and ignoring the evanescent 
terms, the spmor 4'(a;) is analogous to that of a quantum 
spin Hall dotP^Ull 


B. Luttinger liquid description 

The linearized spectrum, obtained in the previous sec¬ 
tion (Eq. allows to construct a LL theory with the 

Fermi energy ep lying near the center of the gap, i.e. 
IcfI ^ ft- Let us introduce the fermionic field operator 
4/(a:) = {if^{x),ifi{x))’^, whose components are 






E 


cie 




+ C 2 e 


— iki 


^>0 


+ -h C4(-l)”e'^3“’(3=-i) 


Cfcl 


with 


( 0 ), 


C4 = (-l)"e"3 ^C4, 


(17) 


(18) 


and a similar equation for 'ifi{x) with Ci —> di- Here 
Cfcj „ is the fermionic operator annihilating the state 
^fel,,^( 2 ^) = ■ In order to con¬ 

struct the Luttinger theory we introduce the right-mover 
field 


V’fl(a:) = 


1 




E 

fci c 


ik-i 


Cki, 


(19) 


with wavevector extended from — oo to -|-oo.^^^f^ As a 
consequence, in the following all the operators must be 
redefined with respect to the vacuum state - i.e with no 
real electrons - | 0 ) by means of the normal-ordering pro¬ 
cedure^ (denoted by : :). From Eq. (11) one verifies 

that 4ir{x) satisfies the twisted boundary condition 

^fl(x + 2 L) = e*"T'^fl(a:), ( 20 ) 

where 7 is given in Eq. (12b |. As for a conventional LLP, 


OBC allow us to express the spinor field in terms of the 
only ifR{x). We have 


lf^{x) = CilpR^x) + C2lfR{-x) 


+ Cge-«3“'-^«(0) + c4e-*t^e'=3"'(-^)^4,(L), (21) 

and similarly for ifiix) with Ci —>■ di- In terms of f’Rix) 
the non-interacting Hamiltonian becomee^ 


( 0 )/ 


Ho = no 


■ '^r{^) {-idx) 'ffR{x) ■■ dx. 


( 22 ) 


l-L 


Following the standard procedure^'0/j(a:) can be written 
via the bosonization formula 


Substituting these coefficients in Eq. (®, the Q-H na¬ 
ture of the states near the center of the gap is evident: in 


IpRix) = 


V2TfA 


giTi-f (Af-|-i7)gz4>(a;) 


(23) 
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with F the Klein factor and N = Ck^ „ : the 

(normal-ordered) particle number operator. Here A is 
the cut-off length, set as A = L/ttNq, with Nq the total 
number of electrons in the dot, and $(a;) the bosonic field 

^iqx-Aql2f ^ ( 24 ) 

Lq 

Here, q = -KTiqlL {riq a positive integer number) and 
bq are bosonic creation and annihilation operators. 

The non-interacting Hamiltonian Hq can be bosonized 



Ho = Vo 

q>0 


(25) 


Let us now introduce the electron-electron interaction. 
As shown in Appendix in the limit L 3 > 1, the 
interacting Hamiltonian has the following form 


Hr,t = 


V{0) 


f-L 


■ [PRix)pR{x)+pR{x)pR{-x)]: dx, 


(26) 

where V (0) is the zero mode of the Fourier transform of 
a short range two-particle interaction and 


PRix) =: iI)\{x)'iI)r{x):= 


N 

TL 


dx‘^{x) 

27T 


(27) 


is the normal-ordered density of right-moving electrons. 
From Eq. (27) and making use of a Bogoliubov transfor¬ 


mation, the total Hamiltonian H = Hq + Ffint can be 
written in the diagonal forrrP 


H = V ^ qdljdq + 

q>0 


2L 




(28) 


where dj, dq are the new bosonic creation and annihila¬ 
tion operators and v = Vo/g, vri = vojg^ are the ve¬ 
locities of bosonic and zero mode respectively. Here, 
Vo is the Fermi velocity introduced in Eq. (12a) and 


g = [1 -|- E( 0 )/ 7 ruo]“^/^ is the Luttinger parameter de¬ 
scribing the intensity of the electron-electron interaction, 
with g < 1 for repu lsive i nteractions and g = 1 for non¬ 
interacting electron^^ZEH jn terms of the new bosonic 
operators dj and dq, the bosonic field in Eq. (24) be¬ 
comes 


5>0 V 9 


:^cos{qx) — ig sin(( 7 x)]d^ -|-h.c.|. (29) 

To be consistent with the linearization assumptions, 
it should be noted that the electron-electron interaction 
must be smaller than half the width of the band-gap. 
This requirement is satisfied when < 70 ( 0 ;, d) < g < with 


flo(a,d) = 


1 + 


aS 

2^271 


'1-f Vl + 4d2 
1 -f 4(52 


(30) 


III. RESULTS 

Let us now discuss the charge and spin properties of 
the dot. Although we stress that our model is general 
and allows to explore the whole range of (5, we will focus 
on the Q-H regime (d < 1/2), where the most striking 
features occur. As already noted, Q-H states lie within 
the gap. We will consider a dot filling such that the 
Eermi surface lies near the center of the band gap. In 
Sec. we will concentrate on characterizing the Q- 

H states. In Sec. |HI B| we will discuss the total charge 
and spin densities, which involve the whole Eermi sea, 
employing a numerical approach based on the Hubbard 
model. 


A. Collective excitations wavefunction and 
magnetization 

Among the most powerful tools to investigate the prop¬ 
erties of the states at the Eermi surface are the collective 
excitation wavefunction (CEWF) 

Pix)= |(Ao|7/>i(a;)|iVo - 1) , (31) 

cr=t,4 

with |7Vo) the ground state with Nq electrons, and 
the collective excitation magnetization (CEM) S(a;) = 
(£ 3 ,( 0 ;), Ey(a;))^ with 

(32) 


^xix) = Yd 


'Fy{x) = 


p=±l 

E 

p—±l 


{Nq\'^\{x) - p'iI}1{x)\Nq - 1 ) 
ipif\{x)\NQ - 1) 


/,tc 


{No\i>\{x) 


,(33) 


while the 2 ; component vanishes identically for simmetry 
reasons, F,z{x) = 0. The CEWE and x component of 
CEM are even functions w.r.t. the dot center, i.e. (p{x — 
L/2) = ip{x+L/2), T,x{x — L/2) = T,x{x+L/2), while the 
y component of the CEM is odd w.r.t. the dot center, i.e. 
'F,y(x — L/2.) = —T,y{x + L/2), with 0 < x < L/2. One 
can directly verify that the amplitude of S(x) is given by 

E(x) = sJf.I[x) F E2(x) = p[x ). (34) 


Equation (31) generalizes the notion of quasiparticle 


wavefun ction, introduced in the context of circular quan¬ 
tum dot j'^° l ^^ l, to the case of a Luttinger liquid. Such a 
quantity probes the probability density of states near the 
Eermi surface and, in the non-interacting regime g = 1, 
reduces to 


p[x) = Lpo{x) = Iff'Aro(x)/ 


(35) 


where 4 ' 7 Vo(a;) is the spinor wavefunction of the A/)-th 
level. Analogously, Eqns. (32 33) describe the magneti¬ 


zation of states near the Fermi surface. For 5 = 1 one 
has 


Yq,{x) = S^q'"\x) = ('^/qJx))'^S-^^jVo(x), (36) 
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with V G {a;, y}. 

Both (p(x) and S(a;) can be extracted via tunneling ex¬ 
periments involving a magnetized STM tip. In particular, 
one has 


(/9(x) cx ^ (x) Vz^, 

(37) 

p —±1 


E,,(x) (X ^ pL^p\x), 

(38) 

p—±l 



where rp‘^^(a;) is the tunneling rate for injecting electrons 
from a tip with spin component p along the quantization 
axis z/. See Appendix [C] for details. 


The Luttinger liquid theory developed here allows to an¬ 
alytically evaluate the CEWF and CEM for any value 
of (5 also in the presence of interactions. Although such 
evaluation can be performed for any temperature, in the 
following we will focus on the most interesting regime 
ksT TTVojL, where expectation values can be per¬ 
formed in the T ^ 0 limit. Useful analytic expressions 
for 6 <C 1/^can also be obtained using coefficients in 
Eqns. 16). Up to first order in 5 one has {x < L/2) 


ip{x) 


^xix) 


T,y{x) 


Kg{0, 0) — 2Kg{0, x) cos{kFx) + Kg{x, x) 

ttA 

—^ [2Kg(0,x)sm(kFx) — Kg(x,x)sm(2kFx)], (39) 
ttA 

2Kg(0, x) sin(fcFa;) — Kg{x, x) sm{2kFx) 

ttA 

r 

—r {-Kg{0,0) + 2Kg{0,x) cos{kFx) 
ttA 

Kg{x,x)[l — 2cos{2kFx)]} , (40) 

Kg{0, 0) — 2Kg{0, x) cos{kFx) + Kg{x, x) cos{2kFx) 
ttA 
2(5 

—^ [2Kg{0,x)sm{kFx) — Kg{x,x)sm{2kFx)], (41) 
ttA 


one has that Wg{x « 0) oc (ttA/L)^/^®, while Wg{x « 
L/2) (X (7rA/L)*^®+^/®^/"‘. Furthermore, for g = 1 one has 
lUi(a;) = (ttA/L)^/^ Vx. 

In the non-interacting case, the above expressions reduce 
therefore to 

ip[x) = [cosh(K 3 °^x) — cos(fci?x)] 

-I- d[2e“”3sin(/ci7’a^) — sin(2fcFx)] |, (44) 

^x{x) = ^ |2e“'‘3sin(fci?x) — sin(2fcFx) 

r (0) 

— 25\^cos{2kFx) — e ”3 ^ cos{kFx) 

-e-«™3^sinh(4°^x)]}, (45) 

Ey (x) = ^ I cos(2kFx) — 2e“'‘3 cos(fc_Fx) -I- 

— 2(5[sin(2fci?x) — 2e“'^3 | ^ ^ 45 ^ 


which c oin cide, as anticipated, wit h the qu antities Pq{x) 
in Eq. (31) and Eq^^(x) in Eqns. (32 33) evaluated for 
X < L/2 to the first order in <5, as can be directly ver¬ 
ified employing the definition of the dot spinor wave- 
function in Eq. ® and neglecting terms proportional to 


exp 


do) 


(x-L) 


as discussed above. 


In the Q-H regime the CEWF and the CEM exhibit 
markedly enveloped oscillations as shown in Fig. for 
non-interacting electrons. We start considering Nq = 49 
and 5 = 0.044 ^ 1, deep in the Q-H regime. 

The CEWF (panel (a)) displays oscillations with a 
wavevector fcp and Np = 25 peaks (taking into account 
that p{x) is even w.r.t. the dot center). For any A^o one 
has 


Np = 


No/2 

{No + l)/2 


if No is even, 
if No is odd . 


(47) 


where we have neglected terms oc exp 
for X < L/2. Furthermore, 


do) 


(x-L) 


< 1 


idg(xi,X 2 ) = (42) 


with 

Wg{x) 


l_2e-A/^cos(^^) 


I g~27rA/L 



(43) 


and kF = ttNq/L tt^/2L the Fermi momentum for 
A’o electrons. The above expressions consist of non¬ 
oscillating terms oc Kg{ 0 , 0 ), terms oc Kg{ 0 ,x) oscillat¬ 
ing with wavevector kF and terms ex Kg{x,x) oscillating 
with wavevector 2k f- Constant terms and terms oscillat¬ 


ing with wavevector kF are enveloped by exp 



X 


and exp 



respectively. 


To lowest order in A/L, 


This behavior is in contrast with the expectations for 
the standard Friedel oscillations of a one-dimensional 
system with non-degenerate stateJ^^^^S®^, which instead 
would predict oscillations at wavevector 2kF, with No 
peaks. We dub this an anomalous Friedel oscillation. 
The CEWF oscillations are more pronounced near the 
edges, while near the center p{x) is flatter. 

On the other hand, the oscillations of Ey(x) (panels (b, 
c)) display No peaks corresponding to a wavevector 2kF, 
in agreement with the expectations for the standard 
Friedel oscillations, with ^^^(x) and Ej,(x) essentially 
out of phase by f/2 - suggesting a precessing pattern 
of the CEM which strongly resembles that of a helical 
systeirPSHH. However, CEM oscillations near the edges 
are far less regular than those in the center. Indeed, 
near the edges an alternating sequence of high and low 
peaks emerges. To better visualize the spin texture, 
Fig.i shows the spatial pattern of the CEM vector in 
different dot regions. Near the edge (panel (a)), where 
the anomalous Friedel oscillations are stronger, the 
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0 X Ll2 



FIG. 3. Plot of the CEWF (a) and of the x, y components 
of the GEM (b, c) for No = 49 non-interacting electrons in 
the Q-H regime with 5 = 0.044. Parameters here: a = 150, 
P = 1000 and g = 1. 



FIG. 4. Vector plot of 5](a;) (a) near the edge and (b) near 
the center of the dot for the Q-H case <5 = 0.044. Gray shades 
denote the quasi-period of the magnetization pattern. Pa¬ 
rameters as in Fig. 


GEM vector length is strongly modulated and an effec¬ 
tive anomalous oscillation of the magnetization, with 
wavevector emerges. This is a direct consequence of 
the anomalous Friedel oscillation. Indeed, the amplitude 
of S(a;) is given by the charge oscillations in g}{x), as 
shown in Eq. (34) . On the other hand, in the center 
(panel (b)) the magnitude of the GEM is only slightly 
modulated, allowing to observe a more conventional 
precession with wavevector 2kp. Therefore, a mixed 
character of the spin pattern emerges, with different 
periodicities between the edges (kp) and the dot center 
{2kp). 



FIG. 5. Plot of the CEWF (a) and of the x,y components 
of the GEM (b,c) for No = 49 non-interacting electrons in 
the Q-H regime with 5 = 0.177. Parameters here: a = 150, 
P = 4000 and g = 1. 


The anomalous Friedel oscillations and the peculiar spin 
textures are a genuine hallmark of the Q-H regime. They 
are driven by the exponential edge terms (oc exp[— 
for X < L/2 and cx exp[K 3 °^(x — L)] for x > T/2), 
peculiar of states in the band gap, with a typical length 


scale I = For a = 150 and P = 0.044 - see 

Fig. ^ - one has k''^^L « 6.5 and thus £ « 0.15L. In 
such a situation, oscillations with wavevector kp fully 
emerge in the GEWF and the peculiar spin texture 
with different periodicities at the dot edges and center 
develops in the GEM. 

To get a clearer picture of how this may occur, let 
us analyze the approximate expressions for !p(x) and 
S(a;) in Eq. ( 44]^ |. For S <C 1/2, as in the case 
discussed so far, the leading term for the GEWF is 
the sole oscillation with wavevector kp. On the other 
hand, the GEM components display both oscillations at 
kp and 2kp - see Eqns. ( 45|46 ). The mixing between 
the two is controlled by the exponential term which 
damps oscillations with wavevector kp near the dot 
center. This mixing is responsible for the irregular, 
alternating pattern of high and low peaks in Figs. |^b, 
c). Increasing S, two phenomena occur. Firstly, 
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' \ i - f 1 y I X ^ i 

0 xjL 0.08 
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0.46 


xj L 


0.54 


FIG. 6. Vector plot of S(a;) (a) near the edge and (b) near 
the center of the dot for 5 = 0.177. Parameters as in Fig. 


increases and consequently i shrinks, resulting in a 
suppression of the anomalous Friedel oscillations of the 
CEWF in the center, accompanied by a shrinkage of the 
region where the mixing between oscillations at and 
2kp occur in the GEM. In addition, corrections to the 
picture valid in the Q-H regime occur, signaled by the 
terms oc i5 in Eqns. ( 44p6 |. As a result, conventional 
Friedel oscillations with wavevector 2kp emerge in the 
CEWF. They are expected to be stroger near the center, 
while anomalous Friedel oscillations still survive in a 
region with length of order £ near the edges. 

Anomalous Friedel oscillations and peculiar GEM 
patterns are expected until £ > L/Nq, the latter being 
the typical wavelength of the conventional Friedel 
oscillations. 

Our predictions are supported by Fig.j^and Fig.|^ which 
show the results for S = 0.177 (with £ « 0.05L > L/Nq). 
The CEWF displays a regular pattern of conventional 
Friedel oscillations in the center and only few anomalous 
Friedel oscillations survive near the edges. Also the 
GEM displays a far more regular behavior, with much 
less evident alternating high-low peaks near the dot 
edges. This also reflects in a more regular precession 
of the vector S(a:) in the latter region, as shown in 
Fig. I^a). Increasing S > 1/2, the dot totally leaves 
the Q-H regime and becomes polarized. Here, the dot 
displays conventional Friedel oscillations while the spin 
begins to become polarized along the x direction (not 
shown) .14^ 


Let us now discuss interaction effects in the Q-H regime. 


with the aid of the approximate expressions in Eqns. (39 


42). For g < 1 one can understand the fate of the anoma¬ 


lous Friedel oscillations inspecting the power-law scaling 
of Kg{xi,X 2 )- Near the dot center, for x « L/2, one has 


^s(0,0)oc^^^ , 

Kg{0,x) oc 
Kg{x, x) oc 


(x)* 


(48a) 

(48b) 

(48c) 


Thus, although near the center both (p{x) and S(x) 
vanish as interactions are increased, anomalous Friedel 



FIG. 7. Plot of the GEWF (a) and of the x, y components of 
the GEM (b,c) for No = 49, 5 — 0.177 and different values 
of the interaction parameter: g = I (blue, solid), g = 0.8 
(red, dashed) and g = 0.6 (green, dotted). Parameters here: 
a = 150, P = 4000. 


oscillations cx Kg{0,x) tend to zero with a faster power 
law than conventional Friedel oscillations oc Kg{x,x). 
On the other hand, near the edges x » 0,L all three 
terms scale with the same power law oc (ttA/L)^/^®. 
Gonsequently, as the interaction strength is increased, 
the region where the anomalous Friedel oscillations can 
be observed shrinks near the dot edges. Our analysis is 
supported by the behavior of the GEWF and GEM for 
different values of g and a fixed value of 6 as depicted in 
Fig. [3 For a given interaction strength the behavior as 
a function of S remains qualitatively similar to that for 
non-interacting electrons already discussed above. 


B. Charge and spin densities 

We now investigate the visibility of the phenomena 
observed previously at the Fermi surface in the Q-H 
regime in the total charge p(x) and spin densities S(x) = 

(^Sxix), Sy{x)^ , which involve the whole Fermi sea. As 

already anticipated in Sec.|llj the ratio between the num¬ 
ber of states in the gap and the total number of states 
in the Fermi sea scales as 6, so that in the Q-H regime 
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with 6 < 1/2 a different approach should be used to 
evaluate th e abo ve quantities. We have employed a Hub¬ 
bard modePsESl the dot has been discretized into Af 
sites and corresponding fermionic operators Cj^a- are in¬ 
troduced. The Hamiltonian iJnub = Htb + H-mt is 

AT-l Af-l 

3=1 i=i 

J ^-1 

i=i i=i 

1 ^ 

+ . (49) 

i=i 

N 

Hint = U ^ (^ 9 ) 

i=i 

Here, t is the hopping amplitude between neighbor¬ 
ing sites and U > 0 the strength of the repulsive 
on-site electron interaction. We are interested into 
p{x) = {p{x)) and S^{x) = {S^{x)), where (...) de¬ 
notes the zero-temperature quantum average valid in 
the low-temperature limit ksT <C ttwo/T. Note that 
Sz{x) = {Sz{x)) = 0 and that p{x) and Sv(x) share the 
same spatial symmetries of ip{x) and Sj^(a;) respectively. 

Let us begin considering Nq = 48 non-interacting 
electrons, which can be treated by means of an exact 
diagonalization. Figure shows the charge and spin 
densities for <5 = 0.044, deep in the Q-H regime. The 
total charge density exhibits oscillations centered around 
the average value Nq/L as expected with Np = 24 peaks, 
consistently with the anomalous Friedel oscillations of 
the CEWF. This number of peaks is also what one would 
naively expect from standard Friedel oscillations con¬ 
sidering that the states below the gap have a parabolic 
spectrum with D = 2. However, the influence of Q-H 
states at the Fermi surface gives rise to a much flatter 
density in the center and more pronounced oscillations 
near the edges. This confirms the picture discussed for 
the CEWF in the previous section. 

Also the spin density S(a:) displays signatures induced 
by the anomalous Friedel oscillations. Both Sx{x) and 
Sy{x) display Nq peaks consistent with a 2 A:f oscillation. 
However, an anomalous alternating pattern of high and 
low peaks is observed in both quantities near the dot 
edges, in strong analogy with the behavior of the CEM 
components. We note in passing that Sx{x) oscillates 
around a non-zero reference level, connected to the 
partial polarization of the dot induced by the external 
magnetic field. Increasing S, the ratio between the states 
in the gap and those below it, increases. However, as 
also discussed above, the kp oscillations induced by 
states in the gap tend to be confined only near the dot 
edges while the dot center becomes dominated by 2kp 
oscillations. As a result, one would expect a mixed 




y (c) components densities (units L~^) for a quantum dot 
with Ao = 48 non-interacting electrons in the Q-H regime 
S = 0.044. The inset in (a) shows the weak density oscillations 
near the dot center. Parameters: L — 8.4/rm, Af = 432, 
t = 4.5 meV, U = 0, a = 150, = 1000. 


character of the charge density, with 2kF oscillations 
in the center and kp oscillations near the dot edges. 
Analogously, one can expect a more regular pattern of 
the spin density oscillations. This is confirmed by the 
numerical results of the Hubbard model shown in Fig. 
for 6 = 0.177. The expected behavior is indeed observed. 

Let us now consider the effects of electron interactions. 
Figure [l0| shows the situation for S = 0.044 and U = t/2. 
Here, calculations are performed via a MPS variational 
algorithnPSl. Oscillations of p{x) with wavevector 2kp 
are now evident in the dot center, while few ones with 
wavevector kp are present in a narrow region near 
the edges. We still attribute this effect to the role of 
the Q-H states near the Fermi surface in shaping the 
electron density. We also notice that the spin density 
exhibits more regular oscillations, as it occurs for the 
CEM when interactions are turned on. Both these facts 
support the idea that even in the interacting regime, 
features of the Q-H states near the Fermi surface can 
be detected in the charge and spin densities of the system. 

We close this section with a brief comparison between 
the results provided by the Luttinger model and those 
of the Hubbard model. A common ground can be found 
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FIG. 9. Same as in Fig. [^but for 5 = 0.177. Here, /I = 4000 
with other parameters as in Fig. 



0 X Ljl 


FIG. 10. Same as in Fig.j^but for U = 112. Other parameters 
as in Fig. 


FIG. 11. (a) Plot of 5p{x) for No = 49, 5 = 0.177 and 

U = t/2 calculated with the Luttinger model (blue, solid) 
and the Hubbard model (red, dashed); (b) same as in (a) 
but for 5Sx{x)-, (c) same as in (a) but for 5Sy{x). Here, we 
have optimized the agreement between the curves by fitting 
the Luttinger parameter g = 0.78. Other parameters as in 
Fig.[Tg 


in the difference between charge and spin densities for 
iVo and iVo - 1 electrons, 5p{x) = pix)\jy^ - p{x)\j.^^_^ 
ajxd SS,^{x) = S'i,(a;)|jY^ — S'i.(a;)|jY^_j^ respectively. These 
quantities can be experimentally probed, e.g. in the shift 
of the chemical potential induced by a charged^ or mag- 
netizecP^ tip. In the non-interacting case, it is immedi¬ 
ate to show that dp{x) = tpoix) and SSi,{x) = Sg‘'\x). 
The results of the Luttinger and Hubbard model in this 
regime coincide within numerical accuracy (not shown). 
In the interacting case the density differences depend es¬ 
sentially on states at the Fermi surface within the limits 
of the Luttinger model developed here. Figure 11 shows 
dp(x) and 6S^{x) for the case U = with Ng = 49 and 
(5 = 0.177. The Luttinger parameter has been fitted here 
to be g « 0.8 by maximizing the overlapping between 
the three calculated quantities. Such value is well within 
the range of validity of our model (see Eq. @). As is 
clear, the agreement between the two models is indeed 
very good. 

































































































































11 


IV. CONCLUSIONS 


We have studied the properties of a quantum dot 
embedded into a spin-orbit quantum wire subject to an 
external magnetic field. Focusing on the regime where 
Q-H states develop in the energy gap, we have built 
an analytical Luttinger theory with open boundaries. 
Both charge and spin properties of the Q-H states have 
been analyzed, focusing on the collective excitation 
wavefunction and magnetization respectively. In the 
Q-H regime they respectively display anomalous Priedel 
oscillations, characterized by a wavevector kp instead 
of the expected 2kp, and peculiar spin textures in 
which the Q-H states magnetization processes with a 
competition between oscillations at wavevectors kp at 
the edges and 2kp around the center. Such effects are 
due to the presence of dot edge states, the hallmark of 
the Q-H regime, which tend to occupy the entire length 
of the dot in the Q-H regime and for weak interactions. 
When the magnetic held is increased, or interactions get 
stronger, the anomalous Friedel oscillations and peculiar 
spin textures tend to disappear from the dot center, 
while they are more robust near the edges. Signatures of 
these peculiar states can be detected also in the charge 
and spin densities although the effects on the latter 
quantities are much weaker. 

The parameters employed in the paper can be achieved 
in state-of-the-art systems. Indeed, for the case of InAs 
nanowires, the typical values of the parameters that ap¬ 
pear in Eq. ([^ ar^SMMl r] m 2 ■ 10 eV ■ m, g* m 5, 
TO* « 0.023 TOe, where rUe is the electron mass. To obtain 
the values a = 150, f3 = 1000 and /3 = 4000 used in the 
paper one has to impose L « 10 um and, respectively, 
B = 0.07 T, B = 0.28 T. 
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Appendix A: Secular equation 


Since we are working in the limit |e| <C /3, Eq. ([^ can be 
expanded to the first order in e//3 (fca = in^): 


ki 


do) 


Kg « 4°^ ( 1 


i+m| 


with 


= fci(e = 0) = ^ 


(0) / ^ 
4 = «3(e = 0) = ;^ 


Vl + 4(52 + 1 


Vl + 4(52 - 1' 


1/2 


1/2 


and 


m = 


y /1 + 452 ’ 


m = -m- 


(A2a) 

(A2b) 

(A3a) 

(A3b) 

(A4) 


The critical condition to solve Eq. (A11 is the large num¬ 
ber of states in the gap, kp — k- tt/T, see Fig. 
For 5 ^ 1/2 it casts into [i/a » 1, while for (5 3> 1/2 
it become s y/d 1. In both regimes it can be seen 
from Eq. ( |A3b ) that the latter are equivalent to impose 
S> 1. Thus, Eq. becomes 


cotan(/ciL) = 


Q; 2 e -I- 2/32 

a2\//32 _ g .2 


+ 2(5, (A5) 


that is 


kiL = + 1) — arctan -|- 28 

« |(2n+1) - -arctan(2(5). (A6) 


Substituting Eq. (|A2a|) and solving for e, we obtain the 

(A7) 


discretized energy spectrum 
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(2n + 1) — arctan(25) — 


with 


Vo = /3/?7ifcJ°^ 


(AS) 


If we now substitute the spectrum of Eq. (A7) in 
Eq. (A2a), we also get the discretization of the wavevec¬ 
tor: 

TT 1 

^ 1 ," ~ bT(^’2 + 1) ~ y arctan(2(5). (A9) 

ZLj hi 

In the end, the spectrum assumes the simple form 


In this Appendix we outline the solution of the secular 
equation in Eq. (|^ 

a^\/— £2 [cos(fciL) cosh(KgL) — 1] = 

+ 2 / 32 ) sin(fciL) sinh(KgL). (Al) 




^0 



(AIO) 


In a similar way one could obtain the discretization for 
Kg but, since Kg appears only in exponential terms (see 
Eq. (§) and |? 7 g£„| <C 1, we will ignore the corrections 
due to discretization and set Kg^ri ~ in all formulas. 
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Appendix B: Interacting Hamiltonian 

In this Appendix we discuss the form of the interact¬ 
ing Hamiltonian of Eq. (26). We start from a general 
interacting Hamiltonian 


to be oc 1/k^PL while (B7c B7dl are oc e “”3 (with 


Hint — X 


/ {x)'^\y)V {x — y)i> {y)'^ {x)dxdy, 

^ Jo Jo 

(Bl) 

with a short range two-particle interaction V{x — y). 
Following the standard g-ology modeP^, we assign dif¬ 
ferent coupling constants to different processes. Insert¬ 


ing in ^(x) the expressions for in Eq. (211, Hint 

can be rewritten in term of tpnix). Terms that cannot 
be rewritten as product of two pn densities (Umklapp 
processes) are damped by fast oscillating exponentials 
exy{±2ikpx), exp(±4zfci?x) and thus will be ignored. In¬ 
troducing the notation Aij = c*Cj + d* d j and using the 
relations An = I and A 22 = I, Eq. (BI) becomes 


fV _ zjextra 

■Hint - H-, + . 


(B2) 


Here 


/v(o) _ £ 

int 2 


-L 


with 


where 


[giPR{x)pR{x) + g2{5)pR{x)pH{-x)\ dx, 

(B3) 

(B4) 


hiS) = 92 -gifiS), 


m = 


26 


1 -b Vi + 452 


(B5) 


We now selP^l^, 92 = Pi = E(0) and gi = V{2kp) = 
£E(0), with £ ^ and define the parameter 


g = 


V{0) 


TTVo 


- 1/2 


(B6) 


All the terms oc prJc) with c G {0,-L} are included in 

extra 


Here we can identify five classes of terms: 


oc / ^pR{0)pR{±x)dx, 


2n^^ (x—L) 


PR{L)pR{±x)dx, 


OC 


oc 


3 '^''3 




4 ^ 3 "^^ (x — 


( 0 ) , 


'pR{0)pR{L)dx, 


^^PR{L)pR{L)dx. 


(B7a) 

(B7b) 

(B7c) 

(B7d) 

(B7e) 


a G {2,4}). Thus, in the limit L ^ 1 one can approx- 

r(0) 


( 0 ), 


imate Hint ~ H^nt- From Eq. (27), the total Hamiltonian 


H = Ho + H^^l , being quadratic in the bosonic operator 
bq, can be written in the diagonal fornP32ll 


H = viS)^qdldq 

q>0 


t^vn{6) 

2L 




Here the velocities of bosonic and zero mode are 

Vo 


v{5) = 
vn{5) = 


Kjsy 

v{^) 

my 



and 


K{6) = 




2 pvo + 94 ^ - 92(6) 
2pvo +94 + g2{6) 


m) 

2pvq 


(B8) 

(B9a) 

(B9b) 

(BlOa) 


(BlOb) 


We neglect here the weak intrinsic dependence of g on 
6 due to Vo, which is particularly weak for 6 < 1 /2, the 
regime of interest in this paper. Equations (BIO) can be 
rewritten as 


K'(6,9,0 = U~'+ {g-^-lf 


K{6,g,0 = 


g-2_i(g-2_l)J(J) 


(Blla) 

(Bllb) 


where f{S) is defined in Eq. (B5|. For g = 1 one has 
9 = K{6, 1,£) = K'{S, 1,£) = 1 VJ, £. Table |T] shows the 
comparison between g, K{S,g,^) and K'{{d,g,^)) foil^ 
£ = 0.1 and for different values of 6. The typical values 
of 9 chosen here fall within the region of validity of the 
theory we develop (see Eq. (30)). One can immediately 
see that g « K{S,g,^) « KyS,g,^). This means that 
the contribution of the term oc gi = V{2kp) is negligible. 
Thus we can set K{S,g,^) = K'{d,g,^) = g in all the 
above relations and observe that they turn into the ones 
that follow from the simplified interacting Hamiltonian 
adopted in the main text (see Eq. (26)). 


Appendix C: Tunneling through an STM tip 


Note that, since Pr{x) varies slowly with x, the terms in 
Eqns. (B7a B7b[ B7e| ) are similar and can be estimated 


Spatial oscillations of states near the Fermi surface can 
be probed experimentally by tunneling of electrons in the 
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TABLE I. Comparison between the values of g, K{5, g, and 
K'{5,g,^) for fixed values of <5 and Here ^ = 0.1. 


9 

6 

K{S,g,i) 

K'{5,g,0 


0.05 

0.800046 

0.799980 

0.8 

0.2 

0.800684 

0.799700 


0.5 

0.803169 

0.798620 


2 

0.811292 

0.795181 


0.05 

0.600090 

0.599915 

0.6 

0.2 

0.601345 

0.598739 


0.5 

0.606224 

0.594252 


2 

0.622140 

0.580502 


linear t ranspo rt regim^^ via a (possibly, magnetized) 
STM tifPSE^, Applying a suitable bias to the tip it is 
for instance possible to inject electrons with a given spin 
direction (which needs not to coincide with the quanti¬ 
zation axis) at a specific location x of the dot. Electrons 
then tunnel through the barriers at the edges and flow to 
drain contacts. In the tunneling limit, when the tunnel¬ 
ing through the tip is the slowest process, it can be shown 
that the linear conductance is essential ly det ermined by 
the tunneling rate through the tip onljEHIlIl From now 
on, we will therefore focus on this rate solely. The tun¬ 
neling coupling between the tip and the dot is described 
by the Hamiltonian 




k,(7 




(Cl) 


where r is the tunneling amplitude, cr =t, j, the spin di¬ 
rection (referred to the z axis), k the wavevector in the 
tip and v the quantization axis of the spin in the tip. 
Furthermore, are operators for electrons in the tip, 
represented on the eigenbasis of (Jz- They are connected 


to the operators on the natural basis of spin eigenstates 
along the direction v., with p = + (p = —) for elec¬ 
trons with spin parallel (antiparallel) to the v axis, by 



1 

+ 

(C2) 

71 ' 

[dtX + d^C) ’ 

(C3) 



(C4) 

71 

[dkX - * 4 '^-) • 

(C5) 


Consider an electron with spin orientation p along the 
axis V tunneling into a dot with Nq electrons at position 
X. The tunneling rate for such a process is given 


r^")(x) = [1 - f{/\ED + AAt)] 

(iVo + i|6[;;(x)|iVo)' 


(C6) 


where 




1 

71 

1 

71 

i 

71 


-Iplix) --Iplix) , 
tplix) + iplix) , 

tp^ix) + itplix) , 

-iplix) - i^lix) . 


(C7) 

(C8) 

(C9) 

(CIO) 


Furthermore, is the density of states of electrons in 
the tip with spin orientation p along the v direction, f{E) 
is a Fermi function and AEjj x sxe the energy differences 
between final and initial dot and tip states respectively. 
By looking up the definitions in Eqns. ( C7|C10 ) it is im¬ 
mediately clear that Eqns. (37|38) hold. 
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